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Abstract 

^' 

Ph| We consider deconfined matter in SU{N) gauge theory as an ideal gas of transversely 

p i' polarized quasi-particle modes having a temperature-dependent mass m{T). Just above 

the transition temperature, the mass is assumed to be determined by the critical behavior 
of the energy density and the screening length in the medium. At high temperature, it 
^ \ becomes proportional to T as the only remaining scale. The resulting (trace anomaly 

\ based) interaction measure A = (e — 3P)/T'^ and energy density are found to agree well 

■ with finite temperature SU{3) lattice calculations. 
(N \ 

1 Introduction 

The quark-gluon plasma in the region Tc < T < 5Tc presents a particularly challenging 
topic of investigation to strong interaction thermodynamics. The most suitable tool for 
^ I such studies is the expectation value of the trace of the energy-momentum tensor, (6^) = 

■ e — 3 P, which measures the deviation from conformal behavior and thus identifies the 
interaction still present in the medium. The aim of the present work is to study the 
temperature behavior of this measure and try to identify the underlying physics which 
causes it. The only ab initio calculations in the range of temperatures of interest here are 
obtained through finite temperature lattice QCD. In particular, pure SU{3) gauge theory 
has been studied extensively, and through finite size scaling techniques the behavior is 
given in the continuum limit [Tj . This case will therefore form the main basis of our study. 

In Fig. [1^ we show the temperature dependence of the energy and the pressure, divided by 
T^, for SU(3) gauge theory, and in Fig. [TJo that of the dimensionless interaction measure, 

MT, ^ |g ^ i^. (1) 

For all quantities, the extrapolation to the continuum limit is shown [T]. In the region 
around and just above Tc, the energy density rises much more rapidly than the pres- 
sure, leading to the observed rapid increase of A(T). Since asymptotically e(T)/T^ and 
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3 P(T)/T^ converge to their common Stefan-Boltzmann value (Svr^/lS), there must be 
some temperature Tp at which the growth rates change roles, with the pressure now in- 
creasing more rapidly. This leads to the peak observed for A(Tp ~ 1.05 T^), followed by a 
somewhat slower decrease. The transition itsself is of first order [2l|3], as expected for a 
theory belonging to the Z3 universality class P]. 
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Figure 1: Temperature dependence (left) of the energy density and the pressure, and 
(right) of the interaction measure, for SU{3) gauge theory p^. 



Two further features of the interaction have become clear more recently. As shown in Fig. 
[2] (left), the decrease of A in the region under consideration here is in good approximation 
given by T~^, so that T^A(T) becomes approximately constant very soon above T^, and 
up to about 5Tc [5J. Moreover, it is seen that A(T) in different SU{N) theories scales very 
well with the number of gluonic degrees of freedom [6]; in other words, A(T)/(A^^ — 1) 
becomes a universal curve, as seen in Figure [2] (right). 



rj 3- 



b2- 



SU(3) - lattice 



Q 1- 



2 3 
T/Tc 



0.4r 

0.3- 
0.2- 
0.1- 

0^ 



SU(3) 
SU(4) 
SU(6) 



1.5 



2.5 
T/Tc 



3.5 



Figure 2: Temperature dependence (left) of (T/Tc)^ A(T) for SU{3) gauge theory [T] and 
(right) of the scaled interaction measure A{T)/{N^ - 1) for SU{N) gauge theories with 
N = 3,4,6. 



In the following, we will first check to what extent any of the observed behavior can be 
accounted for by conventional or modified perturbation theory. Next we turn to the non- 
perturbative approach obtained through the relation of the interaction measure with the 
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gluon condensate and the corresponding bag pressure. Finally, in the main part of our 
work, we shall then show that a quasi-particle approach, based on massive gluonic modes, 
can indeed provide an excellent description of all the features observed in numerical finite 
temperature lattice studies of SU{3) gauge theory . 



2 Weak Coupling Approaches 

In perturbation theory, the interaction measure for pure SU{N) gauge theory is to leading 
order given by [HI?] 

Apert = ~ ^ N^9\T), (2) 
P 288 127r2 ^ ^ ^' ^ ^ 

with 

The perturbative interaction measure thus does show the observed scaling in A^^ — 1 just 
mentioned, assuming Ng"^ is kept constant. In eq. ([3]), A^- defines the lattice scale, which in 
the mentioned SU{3) lattice studies was found to be determined by Tc/A^ = 7.16±0.25. 
In this case, we thus obtain 

_ _U_ 4 _ 47r2 1 1.2 

~ 48^^ ~ '33 {ln[7.16(T/T,)]}2 " {ln[7.16(T/T,)]}2 ' 

At T/Tc = 3, we have Apert — 0.13, which is still about a factor 3 below the (continuum 
extrapolated) lattice result Aiat — 0.4. Hence at this temperature, leading order pertur- 
bation theory cannot yet reproduce the plasma interaction. Nevertheless, we have here 
as = g'^/4:n ~ 0.19 for the strong coupling a^, so that in principle perturbation theory 
seems to be applicable, and we could expect that at somewhat higher temperatures, above 
T/Tc ~ 5 — 10, the perturbative form might account for the lattice result. 

The evaluation of higher order perturbative terms has, however, shown that this is not 
the case. Infrared divergences in finite temperature field theory limit calculations to a 
finite order in the coupling g [8]; for the pressure, the highest perturbatively calculable 
order is g^, and calculations have now been extended to this order |9]. In Fig. [3l we 
show the result of expansions in different order g^ for the pressure in SU{3) gauge theory, 
normalized to the Stefan-Boltzmann limit [lOj. It is seen that in the temperature region 
of interest here, T < 10 Tc, the different orders lead to strong fluctuations; the final form, 
up to and including 0{g^), still considerably undershoots the lattice results. 

Moreover, for an understanding of the interaction effects, a comparison of lattice and 
perturbation theory results for the pressure is in fact quite misleading, since the major 
part of the pressure is given by the ideal gas component. To concentrate on just the 
interaction effects, we return to the interaction measure A(T), and here perturbation 
theory breaks down completely. The next-to-leading order (NLO) form for SU{3) gauge 
theory, 

11 



1 4 
D TT 



(5) 



3 




Figure 3: Perturbative expansions of the pressure in SU{3) gauge theory [9llT0]. compared 
to the finite temperature lattice results [I]. 

becomes positive only for g^>0.27, which with the two-loop form of the coupling, 

9-' = ^ HT/At) + ln[2 HT/At)], (6) 

requires inconceivably high temperatures, above 10^ T^. We conclude that the interaction 
of the plasma in the region of interest here, up to some 10 T^, must definitely require some 
non-perturbative features. 

This situation has triggered numerous efforts to modify the perturbative approach to 
include such features. In one approach [10], the 0{g^) term in the pressure is evaluated 
by a non-perturbative scale determination, using lattice results for magnetic screening. 
This leads to a systematic effective field theory, for which in principle all orders can be 
calculated. Another possibility is given by including sums over certain graph classes, thus 
effectively shifting the point about which the perturbation expansion is performed jTT] . 
In particular, this is studied for the terms dominating at high temperature (hard thermal 
loops, HTL) [T2l[T3] and leads to an improved convergence of the perturbation series of 
the pressure. Both approaches have in common 

• a rather good description of the pressure for temperatures above 3 -5 T^, but 

• the range below about 3 Tc is still not well accounted for. 

• Moreover, the strong order-by-order fluctuations for the interaction measure cause 
some doubt that the last order considered is really close to a "flnal" result. 

To illustrate this, we show in Fig. H] (left) the behavior obtained with the help of a partially 
non-perturbative 0{g^) term [10], and in Fig. H] (right) corresponding results from modifled 
HTL calculations [12], in both cases compared to the form obtained in SU{3) lattice QCD. 
The latter show for e — 3P a decrease as 1/T^, so that T^A(T) becomes approximately 
constant above T^. We see in Fig. H] (right) that in leading (LO) and next-to-leading order 
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(NLO) the breakdown of perturbation theory persists also in a HTL approach, and even 
the inclusion of a partially non-perturbative NNLO contribution cannot reproduce the 
lattice result, neither in size nor in functional form. Such a conclusion had been reached 
before, see e.g. [H]. Recent studies [l5j have shown that in the case of full QCD with light 
quark flavors, the discrepancy between lattice data and weak-coupling results is reduced, 
with quite good agreement down to about 2 - 3 T^, however, neither the approximate 
behavior of A(T) in the range from Tc to about 5 T^, nor the sudden drop in the critical 
region can be thus obtained. 




Figure 4: (T/Tc)^A(r) as predicted in systematic effective field theory (left, [10]) and in 
HTL resummed perturbation theory (right, compared to the continuum extrapola- 

tion of lattice studies [1] . 

In general, the breakdown observed in any perturbative treatment as we enter the tran- 
sition region is of course not surprising. Critical or even pseudo-critical behavior with an 
increasing correlation range is simply not a perturbative phenomenon. We therefore have 
to find a non-perturbative approach to address the behavior of the plasma in this region. 

3 Bag pressure and Gluon Condensate 

One of the earliest attempts to account for the essential non-perturbative features is 
provided by the bag model [THIIIT]. Here one implements confinement in an ideal gas 
picture by introducing a bag pressure, measuring the "level difference" between the phys- 
ical vacuum and the ground state in the colored world of QCD. For the corresponding 
thermodynamics, this means that to the ideal gas partition function, Zo{T, V), a bag term 
is added, 

T\nZBiT,V) = T\nZo{T,V) - BV (7) 

which in principle can be determined from a bag model description of hadron spectroscopy. 
The bag pressure simulates a form of interaction [18] , as best seen by the resulting inter- 
action measure, 

^T) = — . (8) 
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We want to check here to what extent this is a viable description of the QGP interaction 
in the region above T^, assuming B to remain temperature-independent. 

The thermal expectation value of the trace of the energy-momentum tensor, {00 = e—3P, 
is related to the gluon condensate, i.e., to the expectation value of gluon term in the QCD 
Lagrangian [T9] . 

where GJJ^ = gF^^ is given by the gluon field of color a in the QCD Lagrangian. The last 
line of eq. (jH]) is obtained using the leading order perturbative beta function, 

/?(») = i^/ + 0(/). (10) 

The value of (G^) = Gq at T = has been estimated numerically, with Gq = 0.012 ±0.006 
GeV^ as "canonical" value [20]. In both analytical [19] and in lattice studies [1], e — 3P 
is normalized to zero at T = 0, so that 

(ep = e-3P = Gg-G|, (11) 

where G^ is the temperature-dependent gluon condensate. In the temperature range 
below Tc, we expect = Gq, so that e — 3P = 0. If the gluon condensate melts above 
Tc, the interaction measure becomes 

so that B = Gq/4 relates bag pressure and gluon condensate. The value for the latter 
given above leads to a bag pressure B^/^ ~ 230 ±30 MeV, which is in reasonable agreement 
with that obtained from hadron spectroscopy as given by the bag model. 

The color summation in eq. runs over the A^^ — 1 gluonic color degrees of freedom, so 
that we can write 



G'' = ^^(^;^^r) = ^^(^^ - m.^Fn, (13) 



where {F^i,F^^) denotes the gluon field contribution per color degree of freedom. The 
scaling of the interaction measure in A^^ — 1 observed for different SU{N) theories is thus 
also in accord with the bag model dependence, keeping g'^N constant. 

If one assumes that G^ = for T > the interaction measure becomes 

A(T) = 91 = 91 (Is.)' ^ , (14) 

using Tc ~ 0.27 GeV for the SU{?>) deconfinement temperature. The lattice data are 
found to decrease much slower and are, as mentioned, in accord with a 1/T^ dependence. 



^It is known that such an assumption is in general not tenable |21j : a description in terms of a 
temperature-independent bag constant must therefore fail eventually. We want to determine here if it 
makes sense anywhere. 
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We therefore compare in Fig. [5] the results for T^A(T) given by the lattice and the bag 
model forms. The bag model naturally cannot account for the structure immediately 
around (the rise to the peak of A(T)), but it also fails in the temperature region above 
Tc. Combining the bag model with some form of weak-coupling appproach can somewhat 
improve the latter, but it can never reproduce the behavior in the critical region. This 
remains true also in various other, conceptually interesting attempts to modify the power 
of the T-dependence of A(T) [22H21]. 
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Figure 5: The temperature variation of A(T)(T/Tc)^ obtained from the bag pressure, 
compared to the corresponding lattice data [1]. 



4 The Quasi-Particle Approach 

There thus remains the task to find a non-perturbative approach which takes into account 
the critical features arising in the temperature region in the range above Tc, as they were 
obtained in lattice studies. We stay in pure SU{3) gauge theory, where, as mentioned, 
the extrapolation to the continuum limit is known [1]. The basis for our considerations 
here is the study of an ideal gas of constituents ( "quasi-particles" ) having dynamically or 
thermally generated masses [251427] . The behaviour of an ideal gas of such massive gluon 
modes provides automatically the observed scaling and also leads to other features in 
accord with the functional behaviour found in SU{N) gauge theories. 

Interpretations of lattice QCD results in terms of a quasi-particle picture have been given 
in many versions [251429] . In our approach, as in [25l[27], we shall include all interaction 
effects in a dynamically generated mass m(T); most other studies maintain in addition 
a temperature-dependent bag constant. Instead, we want to relate the behavior near Tc 
to the critical behavior of the correlation length, causing the effective mass to increase as 
T ^ Tc from above [25] . 

The partition function of an ideal gas of constituents of mass m{T) is in the Boltzmann 
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limit for SU{N) gauge theory given by 



InZ(T) = 2 



{N^ - 1)V 



dp p2 exp(--7p2 + ^2) = 2^ ^-^ K^im/T), 



(15) 

where -ft'i(x) denotes the Hankel function of imaginary argument. The resulting pressure 
becomes 



,(iV2-l)T 



P{T) = T 
1 



dV 



T 



dp p" exp(-- V'p2 + ^2 = 2^ i^2(W^) 



(16) 



while the energy density is found to be 

T2 (d\TlZ{T)\ 



e(T) 



dT I. 



27r2 



y dp p^ exp(— — + m^) I y^p^ + m? — T 



m 



,(Ar2_ 1)^2^2 



37^2 (m/T) + 



T 



/ dm\ 
yjp'^ + m2 y (iT y 

i^i(m/T)| 



(17) 



In these expressions, we have maintained two spin degrees of freedom for the "massive" 
gluons; we return to this point shortly. Both energy density and pressure thus fall below 
the Stefan-Boltzmann limit, as is observed in the lattice data shown in Fig. [H The 
resulting interaction measure is given by 



A(r) 



,(iV^ 



27r2T4 



dp p^ exp 



1 



T 

{N^ - l)m2 
27r2T2 



p2 _|_ ^2~) ) Jp2 _|_ ^2 _ 3 2^ _ 2^ 



dm 



^Jp)^ + m2 \ dT 



m 


1 dm\ 


T ~ 





iri(m/T) 



If m is A^- independent, the scaling in A^^ — 1 is evident. Moreover, if the effective mass 
m is linear in T, as in any conformal theory, A(T) vanishes. Given a running coupling, 
with = N g^{T)T^/3, we get 



^^j.^_.^ {N'-l)m' ^fd{g'Ny 



127r2T2 



dT 



(19) 



We note, however, that such a "naive" quasi-particle description with finite masses seems 
to encounter a conceptual problem. Physical constituents of non-vanishing mass should 
have three, rather than two spin degrees of freedom, since a longitudinal polarization is 
excluded only for massless particles. The resulting changes in all thermodynamic quanti- 
ties - e.g., the increase of the ideal gas energy density e/T^ from Svr^/lS to 127r^/15 - are 
definitely in disagreement with the observed high temperature lattice results. A simple 
shift to massive gluons thus cannot satisfactorily explain the interactions apparently still 
present in the high temperature gluon gas. More generally, a gauge invariant theory does 
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not allow massive physical gluons; the mechanism leading to effective thermal masses 
must thus be more subtle. The mentioned modified HTL perturbation theory approach, 
in which each order already includes some aspects of gluon dressing, not only leads to a 
rather rapid convergence of the expansion; here the contribution of longitudinal gluons 
vanishes in the limit (? — )■ 0, so that one also obtains the right number of degrees of free- 
dom for the Stefan-Boltzmann form [13] . Moreover, it has recently been argued [3QJ that 
masssive gluons should in fact be transversely polarized, since two massless gluons cannot 
combine to form a longitudinally polarized massive gluon [21]. It thus seems justified to 
use the thermal mass scenario outlined above to address the temperature behaviour of 
the quark-gluon plasma. 

The form of the effective mass entering in a quasi-particle approach description has been 
an enigma for quite some time. At sufficiently high temperature, T remains as the only 
scale, so that there we expect m ~ T. From perturbation theory one obtains in leading 
order for SU{N) gauge theory a thermal screening mass ~ N g'^{T) T/3, but in 
view of the above mentioned difficulties, it seems best to leave the proportionality open. 
As we approach the critical point, perturbation theory in whatever form ceases to be 
applicable. We now have a medium of strongly interacting gluons, and the range of the 
forces between them becomes larger and larger as we approach the critical point. This 
range is governed by the correlation length, or in other words, by the distance up to which 
a given color charge can "see" other color charges. This distance is the QCD counterpart 
of the Debye screening radius in QED; we write it as rniT) = 1/ fi{T), where /i(T) denotes 
the corresponding screening mass. It corresponds to the shift from l/k"^ to l/(/c^ + /x^) 
experienced by the gluon propagator due to the presence of the medium. 

The perhaps simplest view thus is to consider the mass of the quasi-gluon in the strongly 
coupled region to be the energy contained in a volume Vcot of the size defined by the 
correlation range, 

rricritiT) ~ e(T)Kor(T). (20) 
In the case of a continuous transition, the critical part of the energy density becomes 

ecrit~(t-l)'-", (21) 

where t = T/T^ and a is the critical exponent for the specific heat. The correlation volume 
(for three space dimensions) can be written as 

V,or = ^7r I drr^T{r,T), (22) 

where 

specifies the correlation function F in terms of critical exponents u for the correlation 
length ^(T) and rj as anomalous dimension exponent. Combining these expressions and 
making use of the exponent equality relating a and u, we obtain 

m,„,(T) ~ (t - 1)1—2-^ = (t _ lyd+v-u). (24) 



9 



for a continuous transition in three space dimensions. For SU{2) gauge theory, the critical 
exponents are given by the corresponding exponents of the 3d Ising model z/ ~ 0.63 
and r] ~ 0.04, suggesting 

m,ruiT)r^it-l)-'-'\ (25) 

This form is correct only in very near the critical point t = 1; for large temperatures, 
^(t) ~ t, so that the overall form expected for the mass of the quasi-gluon becomes 



m{t) ~ a{t - ly^-'^^ + bt 



(26) 



where a and b are constants. The resulting behavior is illustrated in Fig. [6] (left). It 
would certainly be of interest to check this form directly through calculations in SU{2) 
gauge theory; unfortunately, there does not seem to exist any lattice study providing an 
extrapolation to the continuum, thus eliminating finite lattice size effects. Older studies 
of e(T) and -P(T) in terms of a gluon mass m(T) [25] did in fact lead to the form shown 
in Fig. El 

For SU{3), the transition is of first order [2l[3], so that all quantitites remain finite at 
Tc and an equivalent form cannot be given. Nevertheless, in all cases we have a strong 
increase of both e(T) and A(T) in some range above T^, and so we shall maintain the 
functional dependence flM ) / f25|) with an open exponent c. The resulting quasi-particle 
mass is thus expected to have the form 



with constants a, b, c. 



m{T) 



(t - 1)' 



+ bt, 



(27) 
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bt ^ 
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Figure 6: Left: Expected behavior of the effective quasi-particle mass m(t) in (left) SU{2) 
gauge theory, and (right) as obtained from a fit of lattice data in SU{?>) gauge theory. 



Using this mass, we now determine the parameters a, 6, c by calculating the energy density 
from eq. (17) and A(T) from eq. (|T9ll . The resulting mass and the corresponding param- 
eters are shown in Fig. |6] (right). The fits to energy density and interaction measure are 
given in Fig. Hand are seen to reproduce both quantities very well. We can thus conclude 
that the gluon plasma in SU (3) gauge theory in the temperature region above T^. indeed 
behaves like a medium of quasi-particles with masses generated through non-perturbative 
thermal effects. 
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Figure 7: Interaction measure (left) and energy density (right) for SU{3) gauge theory, 
compared to a quasi-particle description; the fitted constants a and b are given in GeV. 



4.1 The Speed of Sound in the QGP 

In a hadronic resonance gas, one finds [32] that the speed of sound drops to zero at the 
critical point defined by the limit of hadronic matter. It does so because any further energy 
increase goes into making more massive resonances, not into momentum and pressure. On 
the QGP side, in the quasi-particle just discussed, the behavior is very similar. As we 
lower the QGP temperature towards the confinement point, the increase of the quasi- 
particle mass has the same effect. In other words, a temperature increase above lowers 
the mass and thus provides more momentum and pressure, causing an increase in the 
speed of sound. 

The speed of sound, defined as 

c^=(^] =^ (28) 

vanishes at Tc for a continuous transition, because the specific heat Cy(T) diverges there, 
while the entropy density s{T) remains finite. In the ideal gas limit, s{T) ~ 4 cqT^ 
and CviT) ~ 12 T^, so that — )■ 1/3. For the temperature-dependent mass fl27|) . the 
speed of sound can be evaluated numerically, using eqns. f|T6l) and f|T7j) . In Fig. [HI we 
show the resulting behavior obtained in our quasi-particle approach, in comparison to the 
SU(3) lattice results. The two forms agree very well and in fact provide the behavior just 
indicated. 



5 Conclusions 

We have shown that the temperature behavior of the interaction measure defined by the 
trace anomaly of the energy-momentum tensor, A(T) = (e — 3P)/T^, is well described in 
terms of massive gluons with only transverse degrees of freedom. The gluon mass m(T) 
increases sharply as T — > Tc from above, due to the rapid growth of the correlation length 
in the critical region. On the other hand, with increasing temperatures and the approach 
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Figure 8: The speed of sound in the quasi-particle approach, compared to the SU{3) 
lattice result [T]. 

to conformal behavior, m{T) ~ T . The combination of these two effects results in a min- 
imum of m{T) around 1.5 Tc, signalling the transition from critical to hot gluon plasma. 
Even the hot plasma, however, remains strongly interacting; weak-coupling studies do not 
reproduce the observed lattice behavior in the region below about 5 T^. 
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